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O ! Abstract 

' The new result on the QSO absorption lines from the VLT-UVES sample is compared 

O , with the past reports on the time-variability of the fine-structure "constant" derived from 

^ . the Keck/HIRES observation, on the basis of an oscillatory behavior of the scalar field 

' supposed to be responsible for the cosmological acceleration. 

l> 

m ■ 1 Introduction 
> 

\ The many-multiplet method, which played a central role in the search for time-variability of the fine- 
structure "constant" a from the Keck/HIRES data of QSO absorption lines QIE], has been applied 



(N 

, recently to the VLT-UVES sample, with more restricted selection criteria [310]. The wavelengths of 

. the lines have been determined by two different ways, using either laboratory wavelengths of Mg II, I 

f~i ■ and Si II lines with terrestrial isotopic abundances (case 1) or with those of the dominant isotopes (case 

I ' 2) [3]. From case 1, which the authors claim to be more robust, they report no evidence of a changing 

O ' . 

^ ' a, in contrast to the previous Keck/HIRES result. Despite the warning that more careful analysis on 

^ ' the likely systematic errors is necessary before final conclusion is reached |2| , it should be worth trying 

to apply some of the theoretical analyses focusing on different features of the observations accepted 
at face value for the time being, in order to provide hopefully a guide in entangling complications 
involved in the phenomenological analyses. 

There has been a class of theoretical models in which a cosmological scalar field (dilaton or 
quintessence) shows itself through a time-dependent fine-structure constant [5|-|15). In some of them 
the scalar field is expected even to be responsible for the cosmological acceleration |S]-|13|. 

On the cosmological side, the scalar field beyond the linearization regime is supposed generally 
to be trapped to a potential superimposed on a smooth background, like an exponential potential, 
thus causing nearly the constant behavior of the sufficient amount of dark energy, the energy density 
of the scalar field acting as an effective cosmological "constant." Due to this trapping process, the 
scalar field a may show a damped-oscillation-like behavior, as a function of the cosmic time i, and 
hence of the fractional look-back time s = 1 — t/<o with tg the present age of the universe. We may 
also assume that the observed coupling strength is proportional to the scalar field, thus expecting 
y{s) = Aa/a x 10^ — KAa{s), with K a constant. 

Another potentially serious issue comes from the Oklo phenomenon 2ni~EI- The detailed analysis 
of the remnants of the natural reactors supposed to have taken place about 2 billion years ago at Oklo, 
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Gabon, West Africa resulted in either the upper bound |Aq;/q;| ;< 1.0 x 10 ^^7], or the nonzero values 
Aa/a - (0.9±0.1) x lO"'^ or - (0.5±0.1) x 10"'^ for dilferent reasons in Refs. [HI and [THj, respectively. 
It has also been argued that it is unlikely that including the effects of the strong interaction enhances 
the bound more than an order of magnitude |18| . Comparing these estimates with the values at the 
level of 10~^ as expected currently from the QSO observations, we face a difference typically of two 
orders of magnitude or more. An oscillatory behavior might provide a possible reconciliation between 
these two kinds of effects. 

As for the suggested virialization process jl5) . also suggested in this connection, more careful 
analyses seem yet to be applied in view of our unique features of the scalar field with respect to the 
structure of matter coupling and the partially finite-range force [71 [201 ■ 

In view of the uncertainties both on the theoretical and the observational sides, we follow a phe- 
nomenological approach to fit the observed Aa/a by the assumed simple damped oscillation for y{s), 
as was attempted in Ref. jSj applied to the 128 data points from the measurement on the Keck/HIRES 
QSO absorption lines We assumed that the Oklo constraint corresponds approximately to a zero 
of the oscillation at Sokio = 0.142, determining the other three parameters by best fitting the QSO 
data. We then found that the 3-parameter fit is nearly as good as the weighted-mean-fit as far as the 
QSO result is concerned [5]. 

We now improve this fit by incorporating a natural zero at s = 0, arising from the fact that Aa is 
the difference of a from today's value and should vanish at s = by definition. This condition is met 
by offsetting the damped oscillation, as implemented JS] by 



limiting ourselves at the moment to the behavior that allows many oscillations. One easily verifies 
that y(s) vanishes at s = and s ~ Sokio- The relaxation time and the period T are both measured 
in units of t^. 

We add that h k. 2.5, T « 0.22 are expected from our typical cosmological solutions 01^, as 
shown in the beginning of Section 4 of Ref. [S] and Table 2 of Ref. ^S] . These values will be referred 
to as "reference values" in the following, though we do not exclude other solutions. 

Equation permits us to compute 



which may be different from the averaged rate of change; the weighted-mean divided by the average 
s, can be compared directly with the laboratory constraint. 




(1) 



where v/s — vi/si = Wokio/sokio = 2ttT ^ with vi determined by 




(2) 




(3) 
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2 The Keck/HIRES result 



We first apply the fit in terms of Q to tlie more recent result of the 143 data points of Ref. The 
least reduced chi-squared Xrd — 1-015 was obtained for a — 0.020, b — 5.5, and T — 1.352, as plotted 
in Fig. 1. By comparing with their x^d — 1.023 for the simple weighted mean y = —0.573 ± 0.113, or 
the 1-parameter fit in terms of a horizontal straight line, our fit is again nearly as good as theirs, the 
same feature as was shown in Ref. |H]. 



Figure 1: The best fit to the QSO data in Ref. 2. The portion of the curve for s < 0.2 is magnified 
by 100 times. Only the binned data are shown for simplicity of presentation, but the actual fitting 
was made for the entire 143 data points. The dashed line is for the weighted mean, the 1-parameter 
fit in terms of a horizontal straight hne y — —0.573, with Xrd — 1-023, with which our 1.015 is nearly 
comparable. 

Note that if we were to extend the horizontal straight line, as drawn by a dashed line in Fig. 1, 
down to Sokio, we would have largely missed the point, contributing an unacceptably large value to 
the chi-squared because the "error-bar" to be used in the calculation is much smaller than those for 
the QSO data, as mentioned before. On the contrary our curve is designed to pass the zero, which 
should be a good approximation to reality for the same reason. 

More detailed comparison does show a decrease of the total amount of x^, but it remains too 
small to pass either of the Akaike or Bayesian information criterion |22) . We emphasize, however, 
that we have introduced more parameters mainly based on a theoretical ground in connection with 
the cosmological acceleration, not simply for a better fit. 

We find that the preferred value of T is much larger than the reference value, ~ 0.22, as explained 
at the end of the preceding Section. Such a small value of T is entirely outside the confidence region 
of 68%, as shown in Fig. 2. We may define the mass nia- — 27t/T. The values of T = 1.352 and ^ 0.22 
mentioned above correspond to m^/Ho ~ 4.75 and ~ 32, respectively. 

To meet the feature of the "observed" broad distribution, T — 0.87 and 1.94 were chosen in 
Ref. corresponding to nia/HQ = 7.25 and 3.24, respectively. They intended, moreover, to 

respect another constraint from the meteorite dating at s « 0.33 |23j by choosing a rather large 6, if 
interpreted in terms of resulting in overly pronounced \y\ toward the high-s end, thus against the 
observation. Nearly the same tendency is found in the analysis of Ref. JT] with mc/Ho = 0.24 — 0.34, 
or r = 26 — 18. In our fit in Fig. 1, we ignored this meteorite constraint according to our own 
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Figure 2: Confidence volume of 68%, represented by three 2-dimensional contours for b = 4.0, 5.5, 7.0, 
respectively, as shown alongside each contour. The cross for b ~ 5.5 is for the absolutely minimized 
Xrdt whereas the two dots are only for the 2-dimcnsional minima for each b. 

argument On the other hand, T ^ 0.5 has been suggested based on A'^ = 4 supergravity in 

Ref. J though including none of the constraints from Oklo and the meteorite dating. 

We find that {a/a)to ~ —3.8 x lO^^^y^^ turns out to be much smaller than the averaged rate of 
change, ~ 0.6 x lO^^^y^^. 

3 The VLT-UVES result 

Quite different features appear, on the other hand, to emerge in the most recent report from the VLT- 
UVES group[3|0]. The result from 23 data points particularly in case 1 is expressed as a weighted 
mean y — —0.06 ±0.06, with x^d ~ 0.95, which may be interpreted as no evidence for time-dependence 
of a, contrary to the Keck/HIRES result 

We attempt the same type of fit in terms of an offset damped oscillator as in Q . By transforming 
the data originally expressed in terms of redshift, z, into a function of s assuming spatially flat 
Friedmann cosmology with to — 13.78Gy, h — 0.70, and fl^ = 0.7, we found the least Xrd ~ 0.53 
for a = -0.050,6 = 3.1, and T = 0.134, as illustrated in Fig. 3, together with Fig. 4 for the 68% 
confidence region. Surprisingly, this reduced chi-squared is even smaller than 0.95 for the weighted 
mean. In Fig. 3, we recognize several data points fitted particularly closely by the oscillatory curve, 
obviously contributing to decrease the chi-squared. 




Fractional look-back time s 



Figure 3: The best fit to the QSO data for case 1 in Refs. 3,4. The portion of the curve for s < 0.2 
is magnified by 2.5 times. The dashed line is for the weighted mean, the 1-parameter fit in terms of 
a horizontal straight line y — —0.06, with Xrd ~ 0.95, compared with which our 0.53 appears even 
improved. 
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Figure 4: Confidence volume of 68% of the fit in Fig. 3, represented in tlie same manner as in Fig. 2. 
Notice a is now negative. 

The decrease 0.42 in x^d implies the increase of the p- value (in the goodness-of-fit test) by ^ 45% for 
the small degrees of freedom 23-3=20. This also results, unlike with the Keck/HIRES data discussed 
in the preceding section, in the decrease of the Akaike and Bayesian information criteria (given by 

+ 2A: and + k\aN, respectively, with k the number of parameters while N is the number of 
data points PH) by 6.3 and 4.0, respectively. These "improvements" may not be significant enough 
to select the oscillatory fit unambiguously because the ideal statistical conditions are unlikely met due 
to possible unknown error sources. On the other hand, we find that the null result a = yields x?d 
as "large" as 1.10, entirely outside any of the contours in Fig. 4. This also implies that our best-fit 
function departs from the null result by 2.6 standard deviations. 

For these reasons we summarize our analysis by saying that the VLT-UVES data for case 1 allows 
a nonzero oscillating a, despite a null result PIE] that is favored by assuming a uniform time variation, 
or the 1-parameter fit in terms of a horizontal straight line. We even argue that choosing a weighted- 
mean fit is tied with a tacit assumption that the true result is close to the uniform distribution. 
On the other hand, the oscillatory fit should be favored as a better starting dependence if there is 
a theoretical reason to prefer an oscillation, which is motivated by a cosmological acceleration, but 
tends to be averaged out by an assumed flat distribution. 

Furthermore, the favored value of T turns out to be "small," though even somewhat smaller than 
the reference value ~ 0.22 expected typically from the two-scalar model, showing obviously a pattern 
significantly different from the "broad" distribution characterizing the Keck/HIRES result. 

Note also that the fitted curve crosses another zero below Sokio, and that {a/a)ta ~ —0.96 x 
lO^^'^y"^ is larger in size, contrary to the situation in the Keck/HIRES result, than the averaged rate 
of change, <, a/a < 1.2 x 10~^^y~^. We point out that the former estimate is already at the level 
recently reached by the upper bound ^ 2.0 x lO^^^y^^ for laboratory measurement |26| . 

We discovered several other fits with larger than the smallest value 0.527 as shown in Fig. 3, 
still considerably smaller than 0.95 for the weighted mean. Only two examples are mentioned here; 
a = -0.000352, b = 10.0, T = 0.156, Xrd = 0-590 and a = 0.156, b = 1.0, T = 0.253, Xrd = 0.777. The 
latter fit is interesting because the obtained T is close to 0.22, the reference value mentioned above. 
Without a zero below Sokio, unlike in Fig. 3, {d/a)tg is even larger; « -|-2.7 x 10^-^^y~^, which might 
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be on the verge of exclusion by the laboratory constraint |26) . 

We also notice that the VLT-UVES result includes the analysis of some sub-samples, giving still 
smaller chi-squared. We focused on the 12 absorption systems listed in the first line of Table 5 of Ref. 
[3] for "single + double (case 1)," applying the same fit as for the full sample. The set of parameters, 
a = —0.066,6 = 2.55, and T — 0.134, yields the minimized Xrd ~ 0.337 again smaller than 0.552 for 
the weighted mean, y = — 0.077 ± 0.101. The curve of the fit looks nearly the same as in Fig. 3. This 
result seems to demonstrate how robust the presence of an oscillation is. 

Having focused on case 1, claimed to be most probably robust, we now apply the same analysis 
to case 2. The weighted mean gives y — —0.36 ± 0.06 with x^d — 1-03. This case turns out more 
like the Keck/HIRES result; a negative y off the zero by 6 sigmas, and apparently broad distribution 
shown in the damped oscillation fit with a = 0.307, 6 = 0, and T — 0.889, as illustrated in Fig. 5, 
entailing x^d — 0-91 only slightly smaller than that for the weighted mean. Note that 6 = does not 
correspond to the absolute minimum of chi-squared, only at the physical boundary 6 > 0. The fits 
with a > result in much larger chi-squared, close to Xrd ~ 2-52 for the assumed null result. 
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Figure 5: The best fit to case 2 in Ref. 3, representing the minimum Xrd edge of the physically 

allowed range 6 > 0. The portion for s < 0.2 is magnified by 10. The dashed line is for the weighted 
mean with y = —0.36. 

It is interesting to note that the presence of mass m„ does not prevent a{t) from falling off smoothly, 
as seen in Fig. 5.8 of Ref. 0, Fig. 1 of Ref. [S] or Fig. 2 of Ref. This represents another 

mechanism for not disturbing the global structure, somewhat different from the one for a much larger 
mass of a arising from the self-mass which makes the scalar force of finite-range of the macroscopic 
order of magnitude i20i , hence leaving detection of possible WEP violation more remote |7| . 

4 Discussions 

Our analysis based on the damped-oscillation fit as suggested by a possible connection with the 
cosmological acceleration showed that the difference between the Keck/HIRES and VLT-UVES (case 
1) results goes beyond a question merely of the presence or absence of the time- variability of a. It 
appears as if the oscillation of the scalar field shows quite different periods of time dependence. Future 
efforts are expected to separate systematic errors focusing particularly upon a better determination of 
the period, thus testing the theoretical models favoring various oscillatory behaviors .[101-1131. 
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In this respect we also add that there is an exceptional example of non-oscillating behavior in the 
simplified version of the two-scalar model [71 1^ obtained by removing the second scalar field |27| , 
though we then lose the aesthetic advantage of the original model, in which no eternal inflation ensues. 
The scalar field continues to fall overriding maxima of the sine-Gordon potential included in addition 
to the background exponential potential until it is finally trapped to one of the minima. Unlike in the 
two-scalar model, a mini-inflation (persistenly as one of the repeated occurrences thus lessening the 
gravity of the coincidence problem) occurs only with a smooth change of the scalar field as illustrated 
in the lower panel of Fig. 6. Note that JIa nevertheless passes through 0.7 as shown in the upper 
panel, enough to cause an acceleration of the universe. We failed to detect any oscillation of a generic 
to the two-scalar model. This choice, still with the scalar field trapped temporarily, might be favored 
if future observations disprove oscillatory Aa/a and if consistency between QSO and Oklo is achieved 
by other means. As shown in Fig. 7, however, the universe today might happen to be entering the 
final eternal inflation, thus featuring an oscillation, basically as in Refs. nni~CSl- Note also that these 
two different behaviors are discussed in somewhat different context in Ref. jl3j . 




58 59 60 61 

Log t 



Figure 6: A solution in the one-scalar model^^ with the potential V(<j) — e~'*^'^[A + m''(l + cos(Kcr))], 
with C = 1.5, A — 1.65 x 10~^,r7i — 0.42, k = 0.351 in reduced Planckian units, shown only around 
the present epoch, io ~ 13.8Gy or Log to ~ 60.2. The scalar field a varies only smoothly (in the 
lower panel), whereas JIa passes through ^ 0.7 (in the upper panel), hence causing acceleration of the 
universe. 
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